Synchronization in networks of mobile oscillators 
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We present a model of synchronization in networks of autonomous agents where the topology 
changes due to agents motion. We introduce two time scales, one for the topological change and 
another one for local synchronization. If the former scale is much shorter, an approximation that 
averages out the effect of motion is available. Here we show, however, that the time required for 
synchronization achievement increases with respect to that approximation in the opposite case. We 
find that this effect is more important close to the continuum percolation transition point. The 
simulation results are confirmed by means of spectral analysis of the time dependent Laplacian 
matrix. Our results show that the trade-oflF between these two time scales, which have opposite 
effects on synchronization, should be taken into account for the design of mobile device networks. 
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After an initial period of characterizing complex net- 
works in terms of local and global statistical properties 
(e.g., attention turned to the dynamics of their in- 
teracting units [H . A widely studied example of such be- 
havior is synchronization of coupled oscillators arranged 
into complex networks. The interplay between topology 
and dynamics is crucial for synchronization achievement 
(see Ref. Q and references therein). In most studies of 
such systems the network has a fixed structure, but there 
are also many interesting scenarios where the topology 
changes over time in various fields such as power trans- 
mission system consensus problem mobile com- 
munication and functional brain networks 0. 

Within the general framework of time dependent or 
evolving networks, we can identify the particular case 
of networks whose nodes represent physical agents that 
move around but interact with each other only when they 
are close enough. Examples include the coordinated mo- 
tions of robots , vehicles @ , and groups of animals [l^ , 
in which cooperative dynamics emerge. Especially, there 
are many examples where synchronization plays a crucial 
role: chemotaxis 11 1, mobile ad hoc networks [T^], and 
wireless sensor networks (TJ]. Despite the importance 
of this topic, prior research on synchronization in time- 
dependent networks of populations of agents has concen- 
trated so far on two special cases: i) where the network 
topology changes fast [l3 - [l7| . and ii) where the popula- 
tion is dense and arranged in a ring [l^. In the former 
case, the fast-switching approximation (FSA) which av- 
erages out the effect of agent motion is commonly used. 
However, for better understanding of synchronization of 
mobile agents and design of an efficient network, it is very 
important to clarify when and how FSA fails. 

In order to study this point, this letter proposes a gen- 
eral framework in which agents perform random walks in 
a two-dimensional (2D) plane. We consider that each 
agent possesses a mobile wireless device whose state 
is characterized by a phase variable, and the phases 



approach one another through the interaction between 
agents within a certain spatial range. This model is 
well suited for communication problems with short-range 
wireless devices. In this letter we show a general mecha- 
nism of failure of FSA when the time scale of local syn- 
chronization is shorter than the time scale of the topology 
change due to the agent motion. Since we need longer 
synchronization time due to this failure, it is an impor- 
tant factor we should take into account for constructing 
an efficient mobile network. 

Our model consists of N agents moving in a 2D space 
(size L X L) with periodic boundary conditions. Each 
agent moves with velocity v. The angle of the ith agent's 
motion is £,i{tk) G [0, 27r], and it changes randomly at 
discrete time steps tk (tk^i —tk= tm). The evolution of 
the zth agent's position is therefore 



y,(tfc -f At) 



yi{tk) - 



-V COS ^i{tk) At 
vs,'mii{tk)At 



mod L 
mod L, 



(1) 



where At < tm- The motion of the agents is diffusive, 
with a diffusion coefficient of 75 ~ v^tm- 

In this paper the dynamics of the oscillators are based 
on the Kuramoto model [l^, which has been applied to 
technological problems recently [23| . The time evolution 
of the phase of oscillator i is represented as 

ip,{t + Tp) = ip,{t) + n ^ sin{ipj{t) - ipi{t)) (2) 

j.dij<d 



+ (j/i ~ Since mobile de- 



where di 

vices emit signals at discrete time intervals, the individ- 
ual phases are updated at discrete time steps of duration 
Tp. Only devices within a distance d of each other can 
interact to approach their phases. 

When the phase difference is small, Eq. ^ can be well 
approximated by the linearized equation 



N 



(Pi{t + Tp) = ipi{t) - (j'^Lij{t)ipj{t), 



(3) 
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which provides interesting hints about the dynamical be- 



havior. Lij(t) = [ki{t)di 



j{t)] is the time dependent 



Laplacian matrix with Cij{t) = 1 if dij < d (i j) and 
Cij{t) = otherwise. ki(t) represents the number of os- 
cillators that are around i within a range d. 

We first note that the instantaneous coordination of 
agents has statistical properties similar to those of a con- 
tinuum percolation. Simply changing d can be enough for 
the system to enter a different topological (static) config- 
uration. The transition takes place at {N ~ l)7rd^/L^ sa 
4.51 (see [2lj and references therein). In this paper we 
fix iV = 100, L = 200, V = 10, and n; ^ 1, which imply 
dc ~ 24.1. Since agents move and the system is finite, we 
do not observe a sudden transition at d ~ dc- 
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FIG. 1: (Color online) Time evolution of the positions and 
pliases of the agents. The colors are mapped to phases on the 
interval [0, 2n], lines are drawn between two nodes when their 
separation is less than d. Each row corresponds to a different 
set of parameters: rp — l, d — 20 < dc for (a), rp=0.01, 
d = 20 for (b), and rp=0.01, d = 40 > dc for (c). 

Starting from random initial phases and positions, 
we let the system evolve according to Eqs. ([T]) and ([2]) 
[2^ . Figure [T] presents time sequences of four snapshots 
for three different sets of parameters. Synchronization 
emerges through motion and intermittent communication 
between agents, even though a single connected compo- 
nent never forms below dc (Figs.[Ha) and (b)). Above dc 
(Fig. [T{c)), agent motion is not necessary but still helps 
the system to reach its final state. 

In all our simulations, the average phase difference 

{Aip) = ^ ,^^^-1) Ej<fcto ~ fk)'^ decays exponentially 
after an initial transient. We can define then a character- 
istic time T in such a way that (Aip) cx e~*/"^and estimate 
T by fitting the numerical data of {Aip). Additionally, 
nx = T/tp stands for the number of phase updates the 
system needs to reach complete synchronization, and its 
inverse defines the system efficiency. Minimizing ut leads 
to more efficient use of the mobile devices' batteries. 

In the literature related to interacting units in time- 
dependent networks, authors often use the fast-switching 



approximation (FSA) [iM^. FSA assumes that the 
topology changes fast enough, and the entries in the con- 
nectivity matrix are replaced by the probability that two 
units are within the interaction range under completely 
random motion (p = nd'^/L^, for d < L/2). The charac- 
teristic time Tps within FSA is expressed as 

Tfs = -rp/ log[l - a{N - l)p]. (4) 

The effect of agents motion is averaged out, and the pa- 
rameters ti and tm do not appear in Eq. It is im- 
portant to note that it makes sense only when the time 
scale of network variations is much shorter than that of 
the interaction. 

Figure dja) plots nx as a function of d for various val- 
ues of Tp. As intuitively expected, there is a decreasing 
monotonic relationship between tit and d. We have also 
drawn a reference line along Tps /tp, which depends only 
on d. For large Tp FSA is very accurate over a wide range 
of d, since the above mentioned condition for FSA is sat- 
isfied. However, we also identify a region in Fig. [H^a) 
where FSA does not hold for intermediate values of d. 
The size of this region increases as rp decreases. 

It is important to note that these deviations from FSA 
take place close to the continuum percolation transition 
point. This is not a sharp transition; instead, T gradually 
deviates from FSA. To gain a broader understanding, we 
plot T/Tps over the Tp-d plane in Fig. [UJb). For large 
enough rp, FSA is very good irrespective of d. 
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FIG. 2: (Color online) (a) Number of updates required to 
synchronize, nr = T/tp, as a function of d. The solid line 
indicates FSA, Tfs/tp (Eq. (g])). The dotted line is Tsl/tp 
(Eq. Q). (b) The ratio T/Tps in the Tp-d plane. In (a) and 
(b) a = 0.005, and dc represents the continuum percolation 
threshold for an infinite system [SJ]. (c) T and nr versus 
Tp with different values of a for d = 5. Filled and open 
symbols of the same color are based on Eq. and the matrix 
product formulation ((8|, respectively. The colored lines in the 
top panel represent Tps for each case. The black line in the 
bottom panel represents Tp^. 

In order to qualitatively explain the deviations from 
FSA. we consider two characteristic time scales: one 
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for clusters to synchronize and another one for break- 
ing apart. The number of time steps (measured in units 
of Tp) for a cluster to synchronize is, to first order in a, 
Us = l/<j\2{d). Here stands for the smallest nonzero 
eigenvalue of the Laplacian of the cluster. On the other 
hand, the number of steps for a single oscillator to leave a 
cluster of size ^(d) (an increasing function of d if c? < dc) 
is rim = ^'^{d)/v'^TMTp. We can then introduce the ratio 



(5) 



s V^TMTp 

which gives us the dominant time scale, where f{d) = 
£_'^{d)X2{d). Note that the topological parameter d (also 
N and L, if they were changed) appears only in /, while 
the parameters related to agent dynamics appear only in 
the denominator. 

It is clear from Eq. that r] decreases if we increase 
Tp. This fact predicts a transition in the dominant time 
scale as we change Tp which is confirmed in Fig. [21 We 
expect the same transition in f{d) by changing d. Our 
numerical results in Fig. [2] suggest that f{d) is an in- 
creasing function of d well below dc- In the following, we 
analyze three different asymptotic behaviors in detail: 

i) The condition r^-Cl holds for small d and large Tp in 
Fig. [2] In this region, the displacement of agents between 
Tp is large. Thus, the network connectivity changes very 
fast before agents synchronize with their neighbors. Fig- 
ure [IJa) shows the evolution of a system under these con- 
ditions. All nodes in the system approach complete syn- 
chronization at approximately the same rate, and non- 
synchronized nodes may become spatially isolated. We 
call this mechanism global synchronization. FSA is accu- 
rate for systems of this type. 

ii) Starting from the previous case, 77 increases when 
either d is increased or Tp is decreased. For d <^ dc, 
r] > 1 implies that the number of time steps required for 
agent motions to rewire disconnected clusters is larger 
than that required for synchronization inside an isolated 
cluster. Local clusters therefore synchronize very easily 
before the topology changes (Fig. [IJb)). The dynamics 
towards complete synchronization is limited by the mo- 
tion of the agents, and we call this mechanism multiple 
cluster local synchronization. Consecutive positions of 
agents are correlated, and the evolution of the system 
depends heavily on the details of the connectivity pat- 
tern. Since FSA neglects such a correlation, it does not 
properly describe the synchronization dynamics in this 
case. Indeed, as shown in Fig. [SJb), FSA fails by orders 
of magnitude. It is clear from Eq. ([5]) that this region is 
broader if Tp is smaller. 

iii) For d ^ dc and 77 > 1 (implying small rp), the 
whole network is connected {single cluster). In this case, 
agent motion is not necessary for the final synchroniza- 
tion. For a static connected network the characteristic 
time is Tstatic = -Tp/log(l - (TA2), where A2 is the 
second-smallest eigenvalue of the Laplacian matrix Q. 
When the network topology changes, its time average 



is a reasonable upper bound for the characteristic time, 
since it can only be improved by the motion of the agents. 
Figure [^Ja) shows that Tsl is not only an upper bound, 
but also a good approximation for small values of Tp. 
When d is large enough for the system to form a complete 
graph, FSA fits again the numerical result. 

In real mobile networks, the signal interval Tp is one 
of the easiest parameters to control. Thus, it is impor- 
tant to study the dependence of T on Tp with other pa- 
rameters fixed. Figure HJc) describes this dependence. 
Decreasing Tp saturates T, while increasing Tp saturates 
ut. Therefore, there exists an optimal value of Tp in the 
intermediate region that simultaneously achieves rapid 
synchronization and high efficiency. This result demon- 
strates the importance to take into account the deviation 
from FSA, since it causes the saturation of T. 

In order to get some analytical insight, we introduce 
the normal modes of the linear dynamics ([3]) . Note that 
Lij changes with time. Let 9i{t) be the normal modes 
corresponding to an eigenvalue A; at time t which satisfies 
= J2iLiUji{t)9i{t), where Uji{t) is the orthogonal 
matrix with a unit determinant. Multiplying the two 
sides of Eq. (jS]) by the transpose C/J(i + Tp) from the 
left, we get 

N 

Oi{t + rp)^J2 Ohn{t)[l - aXm{t)]0m{t). (7) 



Here, Oi,n{t) = J2i Uj^{t+Tp)Uim{t) is orthogonal. Then 
after an arbitrary number of time steps we get 



ei^{t + nTp) = W 



9=0 
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i„=i 



Ci,+ii,(i 



K{t), (8) 



where Iq denotes the suffix corresponding to an eigen- 
mode at time t+qTp. The product of these matrices sepa- 
rately describes the transformation of the normal modes 
of instantaneous networks by Oi^_^-^i^ and the decay of 
each eigenmode by (1 — crXi^). 

This product of n matrices in Eq. ([8]) can be di- 
agonalized. Let its eigenvalues be (1 — crAi)" with 



= Ai < A2 < 



< Aat, the limit of A^ exists for 



Tsi 



-Tp/(log(l - crA2)) 



(6) 



71 — !■ 00. The characteristic time can then be written as 
T = — Tp/ log(l — crA2). Figure[2l^c) compares T obtained 
by this method to that directly measured from simula- 
tions. Their agreement is excellent, even for smaller val- 
ues of Tp where FSA docs not hold. Our procedure can 
be generalized to any other evolving network (23| . 

The two non-trivial behaviors obtained in the simu- 
lations can also be distinguished in the matrix product 
formulation. For 77 ^ 1, where FSA holds, the agents 
move a sufficiently long distance during rp, and Lijit) 
and Lij{t + Tp) are regarded as independent. Thus, their 
eigenvalues, Xi{t) and A/(i-|-rp) are uncorrelated. Then 
we can expect n^'^i(l - crA/J « e"<'°8(i-'^^)> if we ne- 
glect the fluctuation, where the bracket represents the 
average over eigenvalues of the Laplacian matrix. Then 
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we get T = —rp/ (log(l ~ ""A)). Since the average eigen- 
value of the Laplacian matrix is the average degree, we 
get (A) = {N — l)p. Expanding the characteristic time 
in powers of cr, we have 

Tp/T = a{N -l)p + Oia^), (9) 

which is equal to Tp/Tps up to the lowest order in a. 

For "q ^ 1 and d dc (multiple clusters), there is more 
than one zero eigenmode of the instantaneous Laplacian 
matrix, and Oi (t) corresponding to a non-zero eigenvalue 
Xi{t) ^ vanishes before the topology changes, imply- 
ing that local synchronization is achieved. Hence, the 
dynamics of the system is governed by the decay of the 
zero-eigenmodes caused by the topological change. Even 
if we increase the number of signals by decreasing Tp , we 
cannot get a further decrease of the synchronization er- 
ror between disconnected clusters. Therefore, we expect 
that T converges to a finite value for Tp — > 0. Since a 
appears only in the non-zero eigenmodes (1 — ctA;) which 
are neglected in our approximation, the converged value 
of T docs not depend on cr cither (Fig. HJc)). 

In summary, we have presented a model of interac- 



tions between moving agents that takes into account two 
different time scales: one related to local synchroniza- 
tion in clusters and the other related to the topology 
change. We have shown that when the second time scale 
is greater than the first, more time is required for the 
system to achieve synchronization than the prediction of 
FSA. This new effect is particularly important, because 
it affects the optimal parameter values in terms of syn- 
chronization time and efficiency. Although our model as- 
sumes purely random motion, it could be easily extended 
to more realistic patterns of motion . Our result sug- 
gests that the interplay between instantaneous topology, 
agent motion, and interaction rules plays an important 
role for the performance of mobile systems such as ad 
hoc networks or sensor networks. 
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